University of Shizuoka 



US-00-03 
May 2000 



Quark and Lepton Mass Matrices 
with a Cyclic Permutation Invariant Form 

Yoshio Koide[| 

Department of Physics, University of Shizuoka 
52-1 Yada, Shizuoka 422-8526, Japan 



As an attempt to give an unified description of quark and lepton mass 
matrices Mf, the following mass matrix form is proposed: the form of the 
mass matrices are invariant under a cyclic permutation (fi—^fz, f<2—*f$, f$—*f\) 
among the fermions /j. The model naturally leads to the maximal mixing 
between and u T , and with an additional ansatz, it leads to the well-satisfied 



Abstract 



relations \V,, 



us 





PACS numbers: 12.15.Ff, 14.60.Pq, 11.30.Hv 
* E-mail: koide@u-shizuoka-ken.ac.jp 

1 



1 Introduction 



The times of the phenomenological study of the quark mass matrices have been over. How- 
ever, as far as the study of the neutrino mass matrix is concerned, the phenomenological 
study still have the meaning. It is still important to investigate the unified description of 
the quark and lepton mass matrices from the phenomenological point of view. 

In such a phenomenological study, the key to the unified description is to find a 
fermion basis on which the quark and lepton mass matrices at the unification energy 
scale take a simpler and beautiful form. Concerning this point, recently, Takasugi and 
his collaborators JI] have discussed a neutrino mass matrix on a very interesting basis, 
which is described by a discrete symmetry Z3. However, in their neutrino mass matrix 
model, the charged lepton mass matrix is given by a diagonal form, so that the model 
cannot give any predictions for charged lepton mass spectrum. And, at present, they have 
not applied their idea to the quark mass matrices. However, their basic idea seems to be 
highly promising for a unified description of quark and lepton mass matrices. 

Stimulated by their idea, in the present paper, we investigate fermion mass matrices 
with the form 

M = aE + bS(6), (1.1) 
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where the matrices E and S(8) have been normalized as TtE 2 =TtS 2 (6)=1. Although 
Takasugi and his collaborators |l| have related the form (1.1) to a Z 3 symmetry, in 
this paper, we require that the mass matrix is invariant under a cyclic permutation 
(fi—>f2,f2—>f3,f3—>fi) where fi are quarks and lepton fields (/ = u,d,v,e). The form 
(1.1) is the most general form which is invariant under the cyclic permutation and which 
is Hermitian. The matrix M given in Eq. (1.1) is diagonalized by the tri-maximal mixing 
matrix 
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where uj = e 2w ^ 3 , as follows 
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Note that any mass spectrum with three families can be described by the form given in 
Eq. (1.1), because the three terms (a-, bcosO-, and &sin#-terms) in Eq. (1.1) are trans- 
formed into the three diagonal matrices, diag(l, 1, 1), diag(2, —1, —1), and diag(0, 1, — 1), 
respectively, which are independent of each other. Therefore, we have the same number of 
parameters as the number of the observable quantities (mass values). Besides, if the mass 
matrices of all fermion sectors are given by the form (1.1), then the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix V will be given by V = V T V^ = 1 (or V = V$). We must assume 
an additional term which will break the cyclic permutation symmetry. Nevertheless, in 
this paper, we would like to emphasis that the mass matrix form (1.1) which is given on 
the cyclic permutation symmetric basis will shed a new light on the unified description of 
quarks and leptons. The purpose of the present paper is not to give a theoretical model for 
the unified description, but to show how we can see many suggestive empirical relations 
for the fermion masses and mixings if we take a basis on which the mass matrices are 
cyclic permutation invariant. 

2 Basic assumptions 

Our basic assumptions for the fermion mass matrices Mf are as follows: (i) The mass 
matrices Mf are given by a bi-linear form 

M f = mlKfK), 

and (ii) the matrices Kf have a cyclic permutation invariant form 

K f = a f E-b f S(e f ), (2.2) 

where E and S(9f) are defined in Eq. (1.2). Here, in contrast to the definition of M given 
in Eq. (1.1), we have changed the sign of the coefficient of S(0) such as the angle Of is 
in the range < Of < n/2 for aj > 0, bf > and m[ < ra{ < m{. In the expression 
(2.2), the substantial parameters are only 6//a/ and Of, because we discuss only the mass 
rations in the present paper. 



(2.1) 



The form of Mf, (2.1), may be interpreted by a generalized seesaw scenario || 




with rriL'xmji'xKf and Mpocl, where F are hypothetical heavy fermions in addition 
to the conventional quarks and leptons /. However, in the present paper, we do not 
discuss the origin of the form of Mf given in Eq. (2.1) and confine ourselves to discuss 
phenomenological aspects of the model. 

Note that if we assume bf/af = 1, we obtain the relation || 

2 

mi + m 2 + m 3 = -(^+7^+ v^) 2 , (2.4) 
o 

which is excellently satisfied by the observed charged lepton masses. 

In Table ||, we give values of bf/af and Of which are evaluated from the observed 
values of m{ and the relations 

d-q m l + m 2 + m 3 2 

R = 3 / . — t ; — t ; — n 2 = 1 + ( b f/ a f) , 2 - 5 



and 

tang/ = v^-y^a _ , 2 g x 

Jffl3 + a/?712 — 2Jm\ 



For comparison, we have given the values bf/af and 0/ for the following three cases: (a) 
the observed charged lepton masses and the running quark masses mf(fi) at ji = rriz, (b) 
the fermion masses m{ (//) at /x = Ax = 2 x 10 16 GeV in a non-SUSY scenario, and (c) 
the fermion masses m{(fi) at fi = Ax = 2 x 10 16 GeV in a SUSY scenario (the quark 
mass values have been quoted from As seen in Table [I], the values bf/af and Of are 
not so sensitive to the renomalization group effects (evolution of the Yukawa coupling 
constants), because those have been determined only from mass ratios. 
We may read the values Of given in Table [1] as 

cos# e = y^, co&B d = \\^j, cos U = \I^, (2.7) 

which give the angle values e = 47.39°, d = 52.24° and U = 57.31°, respectively, and 
give the relation cos 2 e — cos 2 d = cos 2 d — cos 2 U = 1/12. Of course, the relations (2.7) 
may be accidental coincidence, and they do not need to be taken seriously. 
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Note that not only the down-lepton masses, but also the down-quark masses give 
bf/a,f ~ 1, so that the down-quark masses also satisfy the relation given in Eq. (2.4). 
[However, the value of m{/m,2 is sensitive to the deviation of bf/af from bf/af = 1. If we 
take bd/cid = 1, we must accept the prediction of mf with the deviation of 2a.] On the 
other hand, in contrast to the values of b e /a e and bd/ad, the value of b u /a u is considerably 
deviated from b u /a u = 1. For a reference, in Table |l|, we list values of tan _1 (6j/aj), where 
we put a/ and bf as aj = cos0/ and bf = sin0y since only the ratio bf/af is meaningful. 
It is interesting that 0/ show (fi u ~ 9d and 4>d — e and e = 45°. 



3 Neutrino mass matrix 

We assume that neutrino mass matrix M v is generated by the seesaw mechanism[§], i.e., 

M v ~ -M D M^Ml , (3.1) 

where Md and Mr are Dirac and Majorana mass matrices, V^M^vr and VrMrUr, respec- 
tively, where v c R = {vr) c = CV^. Since we consider that in a similar way to the charged 
lepton sector, the Dirac mass matrix Mn is diagonalized by the tri-maximal mixing matrix 
Vt as 

V T M D V} = D D = diag(mf , , ) . (3.2) 
The neutrino mass matrix M u given in Eq. (3.1) is written as 



M' v = V t M v Vt ~ -D d (VtM r V±)- 1 D d . 



(3.3) 



If the Majorana mass matrix M# is simply given by Mr = m^l, then, by using the 
relation 
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we obtain the form 
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When we take (mf ) 2 ~mf mf , we obtain m^ ~ jm^l = (m^l = mf mf /m R and sin 2 2# 2 3 = 
1. Thus, we can obtain a natural explanation of the maximal mixing between z/ M and v T 



which is suggested by the atmospheric neutrino data J7|. However, the present oversim- 
plified scenario (M R = mjjl) cannot give |Am 2 2 | <C jAm^l (Am-, = (rrt^) 2 — (m") 2 ) and 
sin 2 2^12 7^ which are required for the explanation of the solar neutrino data f|. In 
order to give a realistic numerical result, we must take other terms in M R and into 
consideration. 

For the structure of M R , an alternative scenario is also attractive : the Majorana 
mass tern M' R = VtM r V^ on the basis u' R = VtVr is again given by the cyclic permutation 
invariant form 

M' R = m R K R . (3.6) 

[However, M' R is invariant under the cyclic permutation not of v Ri , but of u R .] Since M' R 
must be symmetric, i.e., {M' R ) T = M' R , the angle parameter 9 must be zero : 



K R = a R E - b R S(0). 
If we take a special case b R /a R = a/2, i.e., 
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The form of M' v in Eq. (3.10) is already discussed in Ref. |§ for the case b v 
the "democratic seesaw" model [Hill 
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and it is known |J that the case b u = — 1/2 gives the maximal mixing between and v T . 

Furthermore, if we consider a special case with = rnf, which arises from Op = 
in Mb = Kr>K D , we obtain a simple neutrino mass matrix of the form 
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(3.14) 



where x = (\/3/2aRmR)mfmf and x/y = mf /mf. This matrix form given in Eq. (3.14) 
has recently been proposed by Ghosal ||TT| on the basis of discrete Z3XZ4 symmetries. 
The mass matrix M' v given in Eq. (3.14) gives the eigenvalues (m^, — m^, — m^), 



m^ = -(^y 2 + 8x 2 + y), -m» 2 = --(yjy 2 + 8x*-y), -m» = -y, (3.15) 



(we have defined as > and m\ > m v 2 > m^), and the mixing matrix 
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where C\i = cos#i2 = y m\ / \m\ + mj) and S12 = sin6 l 12 = y / \m\ + m^) . There- 
fore, for x^$>y (\mf\^>\m^\), we obtain the relations Am^ ~ 2\[2xy and Am^ ~ 2x 2 
together with C12— S12— 1/V2, which leads to a bi-maximal mixing. Thus, we can give 
a reasonable explanation both for the atmospheric and solar neutrino data regarding 
Am?, ~AmL„„~10~ 10 eV 2 and AmL~Am 2 ~1(T 3 eV 2 . 



4 CKM mixing matrix 

If we apply the matrix form given in Eq. (2.2) to the up- and down-quark sectors, we 
obtain the wrong result V = 1. Therefore, let us modify the matrix Kf by adding a small 
term c/P w : 

K f = a f E-bfS(e f ) + CfP u . (4.1) 
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Poj = ^diag(cu,^ 2 , 1) . 



(4.2) 



Then, the form (4.1) is not invariant under the cyclic permutation (f\ — > f 2 , fi — > f3, fs — > 
/i). The term (urf lL fi R +uj 2 J 2L f 2R +~f 3L f?, R ) is transformed into uj 2 (ujJ 1L f 1R +uj 2 J 2L f 2R + 
f ZL fz R ) under the cyclic permutation. However, if the E and S terms are absent in the 
matrix Kf, the form of P^ is invariant under the cyclic permutation, because the common 
phase factor is unphysical. We regard the small P^ term as a "form" invariant term under 
the cyclic permutation in the extended meaning. 

The modified matrix Kf given in Eq. (4.1), is transformed by the tri-maximal mixing 
matrix V T as 



K) = V T K f V} = ±= 
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where \{ are the eigenvalues of the matrix y/3[afE — bfS(9f)] and they are explicitly 
given m Eq. (1.5) with m 8 -> A</-\/3- 

The mixing matrix Ul is obtained by the diagonalization of the Hermitian matrix 



Mf = m K' f (K' f ) 



/ X 2 X 2 c f X lC } \ 

\ 2 c* f \ 2 2 \ 3 C f 
\ Aic/ A 3 c} X 2 J 



+ |c/| 2 l 



(4.4) 



For a small c/ and A3 3> X 2 , 3> Af , the eigenvalues of M'f are given by 



mi~^m Ai , m 2 ~^m (X 2 2 + \c f \ 2 ) , m 3 ~^m (\ 2 3 + 2\c f \ 2 ) , 



and the mixing angles are given by 



(4.5) 



A 2 — A 1 A 2 



(4.6) 



A3 — A l A3 



(4.7) 



and ^13 ~ 0. (If Cf is complex, i.e., argc/ 7^ 0, the elements U\ 2 = sin $12 and U 23 = sin #23 
are replaced by U i2 = sin Q 12 exp(i arg c/) and f/ 2 3 = sin #23 exp(i argc/), respectively.) 
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If we consider that the value of |c/| 2 in the up-quark sector is sufficiently small so that 



6>" ~0, we obtain 
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0.176, 



(4.f 



which is in good agreement with the experimental value |12[ | V^>/V^ s | = (0.0373 ± 
0.0018)/(0.2205 ± 0.0018) = 0.169 ± 0.008. Note that in order to obtain the relation 
(4.8), the assumption (i) given in Eq. (2.1) is essential. 



When we denote as 



M' f = M° f + -\c f \ 2 l 



(4.9) 



as a trial, let us suppose that the lowest one of the three eigenvalues of the mass matrix 
Mj except for the common constant term |c/| 2 l is assigned to zero, i.e., 



In this scenario, we obtain 



detMj = . 



M 2 -a 2 , 



(4.10) 



(4.11) 



together with m[ = |cy| 2 /3. Neglecting the mixings in the up-quark sector, we obtain 



IK. 



M/Ajj ~ \jm d /m s 



0.224, 



(4.12) 
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which are in good agreement with the experiments. However, at present, there is no 
theoretical understanding behind the assumption given in Eq. (4.10). 

The existence of the term with |c/| 2 ~(A{) 2 slightly changes the predicted values 
of m{ from the case of Cf = 0. Especially the values of m 2 is sizably changed as m{ = 
(A2) 2 /3— >m{— [(A2) 2 + |c/| 2 ]/3. However, for the charged leptons, the existence dose not so 
badly spoil the relation given in Eq. (2.4), because |c e | 2 /(A2) 2 ^(A^/A2) 2 ~(m^/m2)~10~ 3 . 
In order to obtain the best fit values of m{ , the values of bf/af and 9f will slightly be 
changed from the values given in Table 1. 



5 Concluding remarks 

In conclusion, we have pointed out that we can see many interesting relations for the 
quark and lepton masses and their mixings if we take a basis on which the form of the 
mass matrices is invariant under the cyclic permutation >f2, f2~^fo, h—>fi)- We have 
assumed that the mass matrices Mf are given by the bi-linear form [given in Eq. (2.1)], 



o 



where the matrix Kf is given by Eq. (4.1). If we put b e /a e = 1 together with c e = 0, we 
obtain the relation given in Eq. (2.4) which is in excellent agreement with experiments. 
We think that the form of the down-lepton mass matrix with bf/aj = 1 is the most 
fundamental one. Although there is at present no explicit theoretical ground for the 
ansatz given in Eq. (4.10), however, if we accept the ansatz, we can obtain the well- 
satisfied relations for \V US \ ~ ^md/m s and \ V c b\ — ^rud/mb- It is also worth while to note 
that the model naturally leads to the maximal mixing between and v T . (This result 
is independent of the assumption given in Eq. (2.1).) If we take the special values of 
the parameters, we can obtain a beautiful neutrino mass matrix (Ghosal's neutrino mass 
matrix) given in Eq. (3.14) which leads to bi- maximal mixing. 
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Table 1: Numerical values of the parameters [and also <pf = tan. -1 (bf /a f)] and Of, 

which are defined by (2.5) and (2.6) in the text. Input values m{(jj) are quoted from 
Ref. [4], where = 2 x 10 16 GeV. The values of fo//a/ and Of have been evaluated only 
from the center values of m{ (/i). 
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